Anomalous diffusion and elastic mean free path in disorder-free multi-walled carbon 

nanotubes 
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We explore the nature of anomalous diffusion of wave packets in disorder-free incommensurate 
multi-walled carbon nanotubes. The spectrum-averaged diffusion exponent is obtained by calcu- 
lating the multifractal dimension of the energy spectrum. Depending on the shell chirality, the 
exponent is found to lie within the range 1/2 < n < 1. For large unit cell mismatch between incom- 
mensurate shells, r\ approaches the value 1/2 for diffusive motion. The energy-dependent quantum 
spreading reveals a complex density-of-states-dependent pattern with ballistic, super-diffusive or 
diffusive character. 

PACS numbers: 73.63.Fg, 73.23.-b, 72.10.Di 



The understanding of charge transport in structurally 
clean systems, with complex and aperiodic long-range 
correlations, has been the subject of intense debate dur- 
ing the past two decades 0. The quantum dynamics in 
most of these systems has been described as an anoma- 
lous quantum diffusion process, related to the multifrac- 
tal nature of the electronic states and spectra |2| ■ These 
unconventional transport mechanisms are incompatible 
with a dominant transport length scale such as the elastic 
mean free path, and the occurrence of a diffusive regime. 
Recently, the discovery of carbon nanotubes has provided 
a whole class of new quasi-one-dimensional systems with 
spectacular effects of topological arrangements of carbon 
atoms on the electronic spectra 0, . The multi- walled 
nanotubes (MWNTs) are intrinsic incommensurate ob- 
jects since, due to registry mismatch between neighbor- 
ing shells, there are very few cases in which the respective 
symmetries of individual shells allow finding a common 
unit cell for the whole object. In most cases, the unit cell 
length (along the nanotube axis) ratio between adjacent 
shells is an irrational number, and the MWNT taken as a 
whole becomes an incommensurate object la. In recent 
years, mechanical and electronic |rj \K. IsL la Hoi ITU 
properties of double- and triple-walled incommensurate 
nanotubes (i-DWNTs and i-TWNTs, respectively) have 
been intensively investigated. Roche et al. first re- 
ported on anomalous diffusion properties of i-DWNTs 
and i-TWNTs by numerical analysis of the wave-packet 
propagation. Signatures of anomalous diffusion were fur- 
ther inferred by the evaluation of the energy spacing dis- 
tribution of energy levels Q, exhibiting Wigner-Dyson, 
Poisson and semi-Poissonian statistics dependingon the 
position of the Fermi level. In a recent work [llj, the 
existence of a finite (energy-dependent) elastic scatter- 
ing rate for electrons in the outer shell of a disorder-free 
i-DWNTs was analytically shown to be a consequence of 
helicity-detcrmincd selection rules for inter-shell tunnel- 



mg HH1. 

Aim of this letter, is to establish the connection be- 
tween spectral and clarify the dynamical properties of 
i-DWNTs, i-TWNTs and in general MWNTs. Our re- 
sults provide a deeper insight on why, experimentally, 
MWNTs typically, but not always 0,03 exhibit normal 
diffusion behavior, with an energy dependent mean free 
path 0] . However, we predict global anomalous diffusive 
behavior in i-DWNTs and in many i-TWNT systems. To 
this extend, first we calculate the spectrum-averaged dif- 
fusion exponent rj, which describes the spread in time of 
an initially localized wave packet by looking at the mean 
square displacement 



X 2 {t) ~t 2r > 
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For this we evaluate the multifractal dimension 1/2 < 
D-i < 1 of the energy spectrum of various incommensu- 
rate, tunnel coupled, i-DWNTs and i-TWNTs, and use 
the relation r\ = D-\ 0,0]. We find a good agreement 
with previous calculations using a wave-packet propaga- 
tion approach 0,0. We find that the exponent r\ strongly 
depends on the chirality and on the number of shells. 
In particular, for fixed number of shells, e.g. two, one 
can find that r\ is closer to the value 1/2, characteris- 
tic of normal diffusion, if at least one shell has a large 
unit cell. Also upon increasing the number of shells the 
diffusive limit 1/2 is more rapidly approached. For ex- 
ample, for the i-TWNT (6, 4)@(17, 0)@(15, 15) is already 
r\ ~ 1/2 such that for such a tube some elastic mean free 
path can be extracted from the wave-packet evolution. In 
contrast, for the (5, 5)@(17, 0)@(15, 15) i-TWNT, with a 
very small unit cell of the (5,5) shell, the diffusion re- 
mains anomalous, and thus no global mean free path can 
be extracted. While the exponent r\ describes spectrum- 
averaged properties, more information about the inter- 
play between density of states and degree of anomaly 
can be extracted from the energy dependent quantities. 
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Thus wave-packet spreading is investigated by solving the 
time-dependent Schrodinger equation and calculating the 
energy-dependent quantum spreading 



£{E,t) 
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v(E)t 



Tt{S(E -H)}v{E)t 



Here S(E — TL) is the spectral measure operator, whose 
trace gives the total density of states n(E), H is the 
system's Hamiltonian, and X(t) is the position operator 
along the tube axis. Finally, v(E) is the group velocity 
at energy E. In some situations £(E, t) becomes time- 
independent such that energy-dependent elastic mean 
free paths can be defined even if the global exponent r\ is 
not 1/2! 

Our starting model to evaluate wave-packet spreading 
is the tight-binding Hamiltonian for a MWNT with M 
shells, where only one pj_-orbital per carbon atom is kept, 
and with zero on-site energies. With a nearest-neighbor 
hopping 70 on each layer n, and hopping (3 between neigh- 
boring layers, it has the form [l8j |: 



A I 



-P £cos(%)e-^|pi">< P ri, 

(n,m) i,j 



where 



and 



h3, 



are sums over nearest shells and 



nearest neighboring atoms, respectively. Moreover, 9ij is 
the angle between the p l ± and pP^ orbitals, and dij de- 
notes their relative distance. The parameters used here 
are: 7 = 2.9 eV, a = 3.34 A, S = 0.45 A @|. An ab-initio 
estimate gives for the intershell coupling /3 ~ 7o/8 [lif . 
Starting from Eq. the spreading of wave packets is 
now evaluated along two different routes. For our first ap- 
proach, we start by observing that the dynamics of wave 
packets in a system is strongly related to the properties of 
the system's energy spectrum. The latter can be gener- 
ally divided into absolutel y co ntinuous, singular contin- 
uous and pure point parts [19j . The wave packets spread 
ballistically if the energy spectrum is absolute continu- 
ous and are localized in systems with purely point-like 
spectra. If the spectrum is singular continuous the wave 
packets anomalously spread [20L I2H . The singular 
continuous spectrum is a multifractal object which can 
be characterized by a set of fractal dimensions D q [23j |. 
Piechon has shown that, at large times t — > 00, the mo- 
ments spread as X q (t) ~ t qDl - q , with D q the fractal di- 
mension of the energy spectrum |l6l| . Then, the spread- 
ing of a wave packet is determined by The motion 
of wave packets in a system with singular continuous en- 
ergy spectrum will be normal diffusive if D-\ = 1/2, or 
anomalous diffusive if 1/2 < D-\ < 1. The energy spec- 
trum of an incommensurate system is usually singular 
continuous, or maintain intrinsic self-similar features [2|. 



Therefore, one may expect some anomalous diffusive be- 
havior of wave packets in i-DWNTs. In order to calculate 
the energy spectrum of an i-DWNT, we approximate the 
irrational ratio of the unit cell lengths of the inner and 
outer shells by one of its convergents, which is the ratio- 
nal number obtained by truncating the continued fraction 
representation of the given irrational number up to a cer- 
tain term 0, 0] That is, we approximate the i-DWNT 
by a commensurate DWNT. As the chosen convergent 
becomes closer to the irrational unit cell ratio, the en- 
ergy spectrum of the commensurate DWNT gets closer 
to that of the i-DWNT under study. We calculate the 
energy spectrum of the commensurate DWNT by direct 
diagonalization of the Hamiltonian in Eq. (J2J . After com- 
puting the energy spectrum, we count the numbers of 
boxes N(l) with length I = AE/2 n needed to cover it. 
Here AE is the range of the energy spectrum and n is a 
positive integer. Thus, the probability density pi(l) is ob- 
tained for each box by calculating the ratio of the number 
of points falling into the box Aj(Z) to the total number of 
points in the data. The set of fractal dimensions D q can 
be defined as 

N(l) 



$>? +1 (0~** D * +1 , asZ 



i=i 

For a simple fractal object, that is, the probability den- 
sity Pi{l) is the same for each box Aj(Z), all general di- 
mensions defined above have the same values Dp. In 
this case, the probability density is Pi{l) = 1/N(l), with 
N(l) ~ 1~ Df . Therefore is D q+1 = D F for any q. Here, 
Dp is the fractal (box-counting) dimension. In general, 
the distribution of points in a fractal object is different 
for different boxes. That is, the probabilities of finding a 
point in the different boxes are different. If q = — 1, one 
has Do = Dp. For q = — 2 the dimension D_i can be 
calculated by 



D-! 



2 f^o 



Pi 



In I 



(3) 



Hence, the dimension Z?_i is extracted from a linear re- 
gression fit to the plot of Ai = — |l n X)i=i i , 7 1 (0 as 
a function of InL As an example, Fig. n snows results 
for the energy spectrum of the i-DWNT (9,0)@(10, 10). 
As discussed in Ref. ^]] , the incommensurability of the 
two unit cells yields a non-vanishing intershell tunnel- 
ing only when helicity-determined selection rules are ful- 
filled, which can occur only if enough sub-bands in each 
shell become populated. In turn, this yields a finite 
lifetime for electrons in one shell due to effective back- 
scattering processes. Here we show that these features 
are also reflected in the energy spectrum which, as shown 
in Fig.^ exhibits a fractal character. We find that within 
the energy range [-12.664 eV, 12.898 eV] the fractal di- 
mension converges to the value n = 0.88, when the in- 



commensurate ratio Ti 



(10, 10)/^(9,0) 



i 



is approxi- 
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mated by the 6th convergent 15/26 obtained by trun- 
cating the continued fraction expression of 1/V3. This 
value is the same as found from numerical wave-packet 
propagation of initially localized wave packets |(J . Here, 
T( n<m ) — a Q y / 3(n 2 + nm + to 2 )/GCD is the length of the 
axial unit cell of shell (n,m), with the carbon-carbon 
bond length ao and GCD being the greatest common di- 
visor of (2n + m) and (m + 2n). 
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FIG. 1: (Color online) Fractal character of the en- 
ergy spectrum of the i-DWNT (9,0)@(10,10) in the range 
[4.04 eV, 4.08 eV]. Inset (a): Plot of -lniV(Z) vs. lnZ as de- 
termined from the energy spectrum. A linear fit to the data 
corresponding to the shortest lengths I yields the fractal di- 
mension D . Inset (b): Plot of Ai = -§ In P, rl (0 vs - 
In I as determined from the energy spectrum. The slope is the 
fractal dimension 

In order to understand the role of chirality, we 
have also investigated the cases of the (6,4)@(17, 0) 
and (5,5)@(17,0) i-DWNTs, differing in the chiral- 
ity of the inner shell. One has T(i7,o)/'? 1 (5 .5) = V3, 
while T , ( 1 7o)/T , (g4) = l/yl9. The difference in the 
two exponents is noticeable: ??(5,5)<&(i7.o) ~ 0.88 and 
?7(6,4)<a(i7,o) ~ 0.82. We attribute the smaller expo- 
nent of the (6,4)@(17,0) i-DWNT compared to that of 
the (5,5)@(17,0) i-DWNT to the larger unit shell mis- 
match. Having in mind the i-DWNT (5, 5)@(17, 0), we 
add an additional outer shell and consider the i-TWNTs 
(5,5)@(17,0)@(15, 15) and (6, 4)@(17, 0)@(15, 15). For 
the first TWNT the diffusion is still anomalous with 
r ?(5,5)@(i7,o)@(i5,i5) ~ 0.88. However, the effect of the 
additional armchair shell (15, 15) is to randomize the en- 
ergy spectrum of the (6, 4)@(17, 0)@(15, 15). To be def- 
inite, already for the approximation 26 : 26 : 15 to the 
ratios 3VT9 : 3 : \/3 we find the value r\ w 0.60, as shown 
m Fig. [3 Notice that, since is defined in the limit 
of box length I — > 0, for the linear fit in Fig. [2] the five 
points corresponding to the smallest lengths I have been 
considered. 

As second calculation route, the energy-dependent and 
time-dependent wave-packet spreading t{E, t) is ana- 
lyzed numerically. The resolution of the time-dependent 
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FIG. 2: (Color online) Plot of Ai = -± In p- x (l) 
vs. In I as extracted from the energy spectrum of the 
(6,4)@(17,0)@(15, 15) i-TWNT. The linear fit to the data 
corresponding to the shortest box lengths / yield the fractal 
dimension Inset: Plot of —\nN(l) vs. In / for the same 

i-TWNT. The slope is the fractal dimension Do- 




E(eV) 



FIG. 3: (Color online) Energy- dependent behavior of the 
spreading £(E,t) for the (9, 0)@(10, 10) i-DWNT, at several 
times t (solid line t = 1000T, dotted t = 5000T, dashed 
t = 6000T, bold t = 7000T, with T = HH/j ). The rescaled 
total density of states is also shown (dash dotted- line) . 

Schrodinger equation is made by expanding the evolu- 
tion operator e~ tHt on a basis of orthogonal polynomi- 
als. This method has been demonstrated to provide an 
efficient real-space computational framework for an or- 
der N algorithm The total length of the MWNTs 
is finite but several tens of /xm-long in order to limit 
boundary effects. For example, the calculations shown 
in Figs. |3 and 0] are performed with a finite length of 
about 24 /im, whereas the unit time step is taken as 
T = 14ft/7o ~ 3.18 fs. Most generally, the anomalous 
quantum spreading is driven by a time dependent power 
law, £{E,t) = i 21 ^)" 1 /v(E), whose asymptotic regime 
gives the conduction regime. When £(E, t) tends to some 
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finite constant at large times, some energy-local elastic 
mean free path £ C \{E) can be extrapolated as £ C \(E) = 
lim t _ >00 £(-B,t). In Fig. and Fig. 0{main frame), the 
energy-dependent length scale £{E,t) is shown at sev- 
eral chosen elapsed times t, starting from wave packets 
that are homogeneously spread in real space with random 
phase on each lattice point. 

For the case of (9, 0)@(10, 10) (Fig. EJ), and energies 
in the charge neutrality point vicinity, a careful analy- 
sis of propagation over several tens of microns reveals 
no deviation from a ballistic-like regime with r](E) = 1. 
Conversely for a small energy window in a high density 
of states region, a very slowly varying length is observed. 
Here the effect of incommensurability becomes more pro- 
nounced due to the enlarged electronic population avail- 
able for scattering. The result of the energy-averaged 
diffusion exponent gives the global exponent r\ = 0.88 
intermediate between the value 1 and 1/2, and in agree- 
ment with prior analysis. 




E{eV) 

FIG. 4: (Color online) Energy-dependent behavior of £(E, t) 
for the (6, 4)@(17, 0)@(15, 15) i-TWNT, at several times t 
(solid line t = 2000T, dashed t = 4000T, dot-dashed t = 
6000T, bold t = 6710T, with T = 14fi/7o). The rescaled to- 
tal density of states is also shown (dotted-line) . Inset: time- 
dependence of the diffusivity £(E,t)v(E) at various selected 
energies. 

In the case of stronger incommensurate systems, such 
as the (6, 4)@(17, 0)@(15, 15) i-TWNT, the tendency to- 
wards energy-dependent anomalous conduction becomes 
even more pronounced. First the region of 2 eV around 
the charge neutrality point remains almost ballistic, 
whereas the rest of the electronic spectrum shows a very 
slower expansion of wave packet in time. The prop- 
agation time runs from 0.4 ps (t = 2000T) to 1.5 ps 
(t = 671 0T), and clearly £(E, t) either shows anomalously 
slow diffusion, or saturates at very short times (inset in 
Fig.0J. Whenever the saturation limit is clearly reached, 
a mean free path £ e i(E) for the whole object can be mean- 



ingfully extracted. £ c \ is found to be around 200 nm to 
400 nm for energies in between 2eV to 3.5 eV 

To conclude, we found that the energy-averaged expo- 
nent rj depends crucially on chirality and number of cou- 
pled shells with global anomalous behavior strongly pro- 
nounced in i-DWNT. When looking to energy-dependent 
properties, we found that even for i-MWNT exhibiting 
rj(E) ~ 1 (ballistic spreading) close to the charge neu- 
trality point, some elastic mean free path can be defined 
in regions of larger density of states where rj(E) ~ 1/2. 
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